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Abstract. The selection effects applicable to the identification of Galactic supernova rem- 
nants (SNRs) at radio wavelengths are discussed. Low surface brightness remnants are miss- 
ing, as are those with small angular sizes (including young but distant SNRs). Several statis- 
tical properties of Galactic SNRs are discussed, including the surface-brightness/diameter 
(E-D) relation. The wide range of intrinsic properties of Galactic remnants with known 
distances, and the observational selection effects, means that the Z-D relation is of limited 
use to derive diameters and hence distances for individual SNRs, or for statistical studies. 
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1. Introduction 

Over two hundred Supernova Remnants 
(SNRs) have been identified in our Galaxy. I 
have produced several versions of a catalogue 
of Galactic SNRs over the last twenty years, 
the m ost recent revised in 2004 January ( Green 
I200I . Here I review some of the statistical 
properties of Galactic remnants based on 
the most recent version of the catalogue. In 
particular I emphasize the importance that the 
selection effects applicable to the identification 
of Galactic SNRs must be appreciated. 



2. The catalogue 

The catalogue of Galactic SNRs exists in two 
formats. Firs t - published as an appendix to 
iGreenl d2004h - the basic parameters (Galactic 
and equatorial coordinates, size, type, radio 
flux density, spectral index, and other names) 
for each remnant. Second, a more detailed ver- 



sion, available on the World- Wide- Web which 
includes the descriptions of each remnant, ad- 
ditional notes and references. The detailed ver- 
sion of the catalogue is available as postscript 
or pdf for downloading and printing, or as 
web pages. The web pages include links to the 
'NASA Astrophysics Data System' for each 
of the nearly one thousand references. Notes 
both on those objects no longer thought to be 
SNRs, and on the many possible and probable 
remnants that have been reported, are also in- 
cluded in the detailed version of the catalogue. 
It should be noted that the catalogue is far from 
homogeneous. It is particularly difficult to be 
uniform in terms of which objects are con- 
sidered as definite remnants, and are included 
in the catalogue, rather than listed as possi- 
ble or probable remnants which require further 
observations to clarify their nature. Since the 
first version of the catalogue was published in 
lGreen1lfT98l the number of identified Galactic 
SNRs has increased considerably, from 145 to 



See http : //www . mrao . cam. ac . uk/surveys/snrs/ 
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23 1 . Much of this increase has been due to the 
availability of large area radio surveys, partic- 
ularly the Effelsberg survey at 2.7 GHz, and 
the MOST survey at 843 MHz, which are dis- 
cussed further below (Section lXTV 

3. Selection effects 

In practice the dominant selection effects ap- 
plicable to identification of Galactic SNRs are 
those that are applicable at radio wavelengths. 
Simplistically, two selection effects apply (e.g. 
lGreenin"99lh . due to the difficulty in identify- 
ing (i) faint remnants and (ii) small angular size 
remnants. 
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Fig. 1. Distribution in surface brightness at 1 GHz 
of 217 Galactic SNRs. The dashed line indicates the 
surface brightness completeness limit discussed in 
SectionEm 



3.1. Surface brightness 

SNRs need to have a high enough surface 
brightness for them to be distinguished from 
the background Galactic emission. This se- 
lection effect is not uniform across the sky, 
both because the Galactic background varies 
with position, and because the sensitivities of 
available wide area surveys covering differ- 
ent portions of the Galactic plane vary. The 
most recent large-scale radio surveys that have 
covered much of the Galactic pl ane are: (i) 
the Effelsberg survey a t 2.7 GHz dReich et alJ 
ll990HFurstet alJl99fi . which covered 358° < 
I < 240° and \b\ < 5°; and (ii) the MOST 
survey at 843 MHz dWhiteoak & A. J. Greer] 
I19961IA. J. Green et alJl999l) 2 . which covered 
245° < I < 355°, but only to \b\ < 1?5. In the 
current catalogue of SNRs, fainter remnants 
are relatively more common in the anti-centre 
and away from b = 0° , where the Galactic 
background is lower (see lGreerl2004l) . 

Since the new SNRs identified from the 
Effelsberg survey were included in t he ver- 
sion o f the SNR catalogue published in lGreenl 
J1991I) . the surface brightnesses of remnants in 
the survey region that have subsequently been 
identified are useful for estimating the com- 
pleteness limit for this survey. Since 1991 an 
additional 24 remnants within 358° < I < 240° 

2 To avoid confusion between authors with the 
same surname, I include initials for Anne J. Green, 
of the University of Sydney, Australia. 



and \b\ < 5° have been included in the cata- 
logue, most in the first quadrant. The surface 
brightnesses of these remnants suggest a com- 
pleteness limit of E sb 10" 20 W m~ 2 Hz ' 1 sr" 1 , 
at 1 G Hz, for the Effelsberg survey (see lGreenl 
2004, for further discussion). And since the 
number of remnants brighter than this value is 
similar in the 1st and 4th quadrants, then a sim- 
ilar limit seems appropriate for the MOST sur- 
vey region. 

Thus, the surface brightness limit for 
completeness of the current catalogue of 
Galactic SNRs is approximately 10~ 20 W 
m _2 Hz _1 sr _1 . Fig. [2 shows a histogram of 
the surface brightnesses of the 217 Galactic 
SNRs, of which 64 are above this nominal 
surface brightness limit. 3 The SNRs with sur- 
face brightnesses below this limit are predomi- 
nantly in regions of the Galaxy where the back- 
ground is low, i.e. in the 2nd and 3rd quadrants, 
and away from b = 0°, as shown in Fig. [2 

Ongoing and future observations will no 
doubt continue to detect more Galactic SNRs, 
although it seems very likely that most of these 
objects will be faint, and hence difficult to 
study in detail. Currently there are several large 
scale radio surveys underway that will cover 
much of the Galactic plane 4 . As is discussed in 

3 For 14 catalogued remnants no reliable radio 
flux density, or only a limit is available. 

4 See: |http : //www. ras . ucal gary . ca/IGPS/ for fur- 
ther information. 
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Fig. 2. Galactic distribution of (top) all Galactic SNR and (bottom) those SNRs with a surface brightness 
at 1 GHz greater than 1CT 20 W irT 2 Hz 1 sr -1 . (Note that the latitude and longitude axes are not to scale.) 




20 40 60 80 100 

Angular size 8 / arcmin 

Fig. 3. Histogram of the angular size of 219 
Galactic SNRs (12 remnants larger than 100 arcmin 
are not included). 



Section l5T2l there is a general trend that fainter 
remnants tend to be larger, and hence on av- 
erage older, than brighter remnants. However, 
because of the wide range of properties of 
Galactic SNRs with known distances, the sur- 
face brightness selection effect applies not just 
to old remnants, but also to young remnants. 



3.2. Angular size 

Current catalogues are also likely to be miss- 
ing small angular size SNRs, as their structure 
is not well resolved by the available Galactic 
plane surveys, and they would not be recog- 
nised as likely SNRs. Fig.[3]is shows the his- 
togram of the angular sizes of known remnants, 
which peaks at around 10 arcmin. The limiting 
angular size varies for the different available 
wide area surveys. As discussed above, the ra- 
dio survey that covers most of the Galactic 
plane is the Effelsberg 2.7-GHz survey, which 
has a resolution of w 4.3 arcmin. So, for this 
survey, any remnants less than about 13 ar- 
cmin in diameter (i.e. 3 beamwidths) are not 
likely to be recognised from their structures 
(although, as discussed below in Section l3~3l 
some searches have been made for small rem- 
nants among compact sources in the Effelsberg 
2.7-GHz survey). The MOST 843-MHz sur- 
vey has a much better resolution, w 0.7 ar- 
cmin, which implies that in the region of the 
Galactic plane covered by this survey only 
remnants smaller than about « 2 arcmin (i.e. 3 
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beamwidths) might be expected to be missed. 
However, although the MOST survey detecte d 
18 new SNRs dWhiteoak & A. J. GreeiJl996l) . 
the smallest new remnant was G345 -7-0-2, 
which is 7 x 5 arcmin 2 in extent, i.e. several 
times larger than the nominal limit of « 2 ar- 
cmin. Thus it is difficult to quote a single an- 
gular size selection limit for current SNR cata- 
logues, although it is clear that it is difficult to 
identify small angular size remnants from ex- 
isting wide area surveys. 

3.3. Missing young but distant SNRs 

The lack of small angular size remnants 
is particularly clear when the remnants of 
known 'historical' Galactic supernovae (see 
Steph enson & Greet] 120021) are considered. 
These remnants are relatively close-by - as is 
expected, since their parent SNe were seen his- 
torically - and therefore sample only a small 
fraction of the Galactic disc. Consequently 
many more similar, but more distant remnants 
are expected in our Galaxy (e.g. lGreenlll985l) . 
but these are not present in the current cata- 
logue. 

Table Ogives the distances, angular sizes, 
flux densities and surface brightnesses at 
1 GHz, for the remnants of known historical 
supernovae from the last thousand years, plus 
Cas A (which although its supernova was not 
seen - so it is not strictly a historical remnant - 
is known to be only about 300 years old). This 
table also lists the parameters of these rem- 
nants when scaled to larger distances of 8.5 and 
17 kpc, i.e. to represent how they would ap- 
pear if they if they were at the other side of the 
Galaxy. 

From Table [2 any young SNRs in the 
Galaxy similar to the known historical rem- 
nants, but in the far half of the Galaxy, would 
generally be expected to have angular sizes 
less than a few arcmin, usually with high 
surface brightness, greater than as 10~ 19 W 
m _2 Hz _1 sr _I (although remnants similar to 
the remnant of the SN of ad 1006 would be 
much fainter). These remnants would also be 
expected to lie close to the Galactic plane, 
with \b\ < 1°. Several dozen other young (i.e. 
less than a thousand year old) SNRs are ex- 



pected in the Galaxy (see lGreenl200 4. for fur- 
ther discussion). However, there are very few 
such remnants in the current Galactic SNR 
catalogue. In fact there are only 3 known 
remnants with angular sizes of 2 arcmin or 
less: Gl-9+0-3, G54-1+0-3 and G337-0-0-1. 
(It should be noted that there are unlikely to 
be other luminous remnants in the Galaxy like 
Cas A and the Crab nebula. Any such rem- 
nants, even on the far side of the Galaxy, would 
have relatively high flux densities, and the na- 
ture of all such sources in the Galactic plane is 
known.) 

Since the missing young but distant rem- 
nants are expected to have angular sizes of 
a few arcmin or less, they will not have 
been resolved sufficiently in large area ra- 
dio surveys. Higher resolution, targeted ob- 
servations are needed to identify any such 
small remnants, and although there have been 
several su ch searches for remnants of this 
type (e.g. iGreen & Gull 11984 iHelfand et al l 
| l985HGreenlll985[ Il989t ISramek et al.lll992t 
Misanovic et al. 2002, see also Saiki aetaLl 
2004), they have had only limited success 
(identifying the small remnants Gl-9+0-3 
and G54- 1+0-3 noted above). (An additional 
candidate small, presumably young SNR is 
G337-2+0-1, which was listed as a possible 
SNR by Whiteoak & A. J. Green. Recently 
ICombi et alJ d2005l) - see also these proceed- 
ings - have also noted that this source is 
associated with an X-ray source listed by 
Sugiz aki et"al] lf2001). which supports the SNR 
identification. 5 ) 

The fact that such missing, small remnants 
are likely to be in complex regions of the plane 
may mean that confusion is a very significant 
problem, and not just at radio wavelengths. 
Further searches for these missing young but 
distant remnants are clearly required. 

4. Some simple SNR statistics 

In the current version of the catalogue, 77% 
of remnants are classed as shell type (or pos- 
sible shell), 12% are composite (or possible 

5 Note that the radio surface brig htness for 
G337-2+0-1 given bv lCombi et"ail <2005l) is too low 
by two orders of magnitude. 
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Table 1. Parameters of known historical SNRs, plus Cas A. 
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0.55 


AD 
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2.3 


8 
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1.1 


1.0 


AD 
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3.2 


7 


1.0 x 10~ 19 


33 
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1.3 


1.2 


AD 
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1.9 


6 


4.4 x 10- 18 


1040 


1.4 


52 


0.7 


13 


AD 


1006 G327-6+14-6 


2.2 


30 


3.2 x 10- 21 


19 


7.7 


1.3 


3.9 


0.31 



composite), and 4% are filled-centre (or possi- 
ble filled centre) remnants. The remaining 7% 
have not yet been observed well enough to be 
sure of their type, or else are objects which 
are conventionally regarded as SNRs although 
they do not fit well into any of the conventional 
types (e.g. CTB80 (=G69-0+27), MSH 17-39 
(=G3577-0-l)). 

There are 14 Galactic SNRs that are either 
not detected at radio wavelengths, or are poorly 
defined by current radio observations, so that 
their flux density at 1 GHz cannot be deter- 
mined with any confidence: i.e. 94% have a 
flux density at 1 GHz included in the catalogue. 
Of the catalogued remnants, 36% are detected 
in X-ray, and 23% in the optical. At both these 
wavelengths, Galactic absorption hampers the 
detection of distant remnants. 

5. Distance dependent SNR statistics 

5.1. Distances to SNRs 

Accurate distances are not available for most 
known SNRs, which is a problem for many 
studies of Galactic SNRs, where it is neces- 
sary to know the distances to remnants (or 
equivalently their physical sizes, since their 
angular sizes are known). The distances that 
are available are obtained from a wide vari- 
ety of methods, each of which is subject to 
their own uncertainties, and some of which 
are subjective. In the next few years further 
distance measurements should become avail- 
able for more SNRs, both because of the var- 
ious ongoing Galactic plane survey, and new 
distance measurement techniques: e.g. from 



Hi co lumn densities (see iFoster & Rou tledge 
or H i absorpti on to polarise d emission 
fro m remnants l|Kothes et al]l2004l) . Currently 
(see [Green 2004), the distances to 47 Galactic 
SNRs are available, i.e. only 20% of known 
Galactic SNRs. The uncertainties in these dis- 
tances are far from uniform. For kinematic dis- 
tances - which are a large majority of the avail- 
able distances - there are always some uncer- 
tainties in deriving distances from observed ve- 
locities, due to deviations from circular mo- 
tion (especially an issue for nearby remnants, 
and for those near I = 0° and 180° where 
the observed velocity does not depend strongly 
on distance) and ambiguities inside the Solar 
Circle. 

5.2. The I-D and L-D Relations 

Many statistical studies of Galactic SNRs have 
relied on the surface-brightness/diameter, or 
'27-D' relation to derive distances for indi- 
vidual SNRs from their observed flux den- 
sities and angular sizes. For remnants with 
known distances (d), and hence known diam- 
eters (D), physically large SNRs are fainter 
(i.e. they have a lower surface brightness) 
than small remnants. Using this correlation be- 
tween 27 and D for remnants with known dis- 
tances, a physical diameter is deduced from the 
distance-independent observed surface bright- 
ness of any remnant. Then a distance to the 
remnant can be deduced from this diameter and 
the observed angular size of the remnant. 

The E-D relation for Galactic SNRs with 
known distances is shown in Fig. |4] This 



6 



Green: Galactic SNRs 



Crab nebula 



X » 



Diameter / pc 

Fig. 4. The surface brightness/diameter (Z-D) re- 
lation for 47 Galactic SNRs with known distances. 
Note that the lower left part of this diagram is likely 
to be seriously affected by selection effects. 



clearly shows a wide range of diameters for 
a given surface brightness, which is a severe 
limitation in the usefulness of the E-D rela- 
tion for deriving the diameters, and hence dis- 
tances, to individual remnants. For a particu- 
lar surface brightness, the diameters of SNRs 
vary by up to about an order of magnitude, or 
conversely, for a particular diameter, the range 
of observed surface brightnesses seen varies by 
more than two orders of magnitude. 

The correlation shown between surface 
brightness and diameter in Fig. |4] is, however, 
largely a consequence of the fact that it is a plot 
of surface-brightness - rather than luminosity - 
against diameter, D. Surface brightness is plot- 
ted, because it is the distance-independent ob- 
servable that is available for (almost) all SNRs, 
including those for which distances are not 
available. For remnants whose distances are 
known, we can instead consider the radio lumi- 
nosity of the remnants. Since E and luminosity, 



Diameter / pc 

Fig. 5. The luminosity/diameter (L-D) relation for 
47 Galactic SNRs with known distances. 



L, depends on the flux density S , angular size 
0, distance d and diameter D, as 



then 



E oc — 

8 2 



E oc 



and L oc S d 



(9d) 2 



or E oc 



D 2 



Thus, much of the correlation shown in the 
E-D relation in Fig. 0] is due to the D 2 bias 
that is inherent when plotting E against D, in- 
stead of L against D. The L-D relation for 
Galactic SNRs with known distances in Fig. [5] 
shows that there is wide range of luminosi- 
ties for SNRs of all diameters. Cas A is the 
most luminous Galactic SNR, but it appears to 
be at the edge of a wide distribution of lumi- 
nosities. The wide range of luminosities is per- 
haps not surprising, given that the remnants are 
produced for a variety of types of supernovae, 
and that they evolve in regions of ISM with a 
range of properties (e.g. density), which may 
well effect the efficiency of the radio emission 
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mechanism at work. For example, some SNRs 
may initially evolve inside a low-density, wind- 
blown cavity, and then collide with the much 
denser regions of the surrounding ISM. 

Furthermore, the full range of intrinsic 
properties of SNRs may be even wider than 
that shown in Figs |4] and |5J as the selection 
effects discussed in Section [3] mean that it is 
difficult to identify small and/or faint SNRs. 

Also, it is not clear that any best-fit E-D 
relation - not withstanding selection effect 
problems - actually represents the evolution- 
ary track of individ ual SNRs (see, for ex ample, 
lBerkhuiisedll98d or iBerezhko & Volkll20o3 
for a recent discussion). The distribution of 
SNRs with known distances is a snapshot in 
time of a population of remnants, and individ- 
ual remnants may evolve in the E-D plane in 
directions quite different from any power law 
fitted to the overall distribution of SNRs (or 
to the upper limit of the distribution). As a 
simple example, consider the situation where 
SNRs have a range of intrinsic luminosities, 
expand with a constant luminosity up to some 
particular diameter - which varies for differ- 
ent SNRs depending on their environment (e.g. 
the surrounding ISM density, which influences 
their expansion speed, which may affect the ef- 
ficiency of radio emission mechanism) - after 
which their radio luminosity fades rapidly. In 
this case, the locus of the upper bound to the 
highest surface brightness remnants for a par- 
ticular diameter is related to where the lumi- 
nosities of different SNRs begin to decrease, 
and does not represent the evolutionary track 
of any individual remnant. 

5.3. A question of regression 

There is a further issue related to any power 
law E oc D" fits to the observed properties 
of samples of SNRs, even if any problems 
with the selection effects are neglected. The 
E-D relation has two particular uses, first, 
the derivation of diameters (and hence dis- 
tances) for Galactic remnants from their ob- 
served surface brightnesses, and second, to 
parameterise the relationship between surface 
brightness and diameter, for comparison with 
models and theories, or between different sam- 



ples of SNRs. For the former, the measured 
value of E is used to predict D, and in this 
case a least squares fit minimizin g the devia- 
tions in log D should be used (see llsobe et alJ 
1990 for a discussion of least squares fitting). 
This, however, has not always been done, and 
instead fits minimizing the deviations in log E 
have been used. The differences between these 
two least square fits can be very large for rem- 
nants much fainter or brighter than the aver- 
age surface of the remnants with known dis- 
tances used to calibrate t he E-D rela tion. For 
example, Case & Bhattacharya ( 1998) derived 
a best fit E-D relation with a slope (n above) 
of -2.38 ± 0.26 using a sample of 36 shell 
SNRs (excluding Cas A) evidently by mini- 
mizing the square of the deviations in \ogE. 
Instead minimizing the deviations in logD - 
which is appropriate if the relation is to be used 
to predict diameters from surface brightnesses 
- leads to a much steeper fit, with a power 
law slope of -3.37 + 0.35. This means that 
Case & Bhattacharya's "best fit" overestimates 
the diameters, and hence also the distances of 
fainter remnants. To quantify this, a faint rem- 
nant with a surface brightness of 3 x 10~ 22 
W m -2 Hz 1 sr~' at 1 GHz, would have its di- 
ameter overestimated by about 40%, which is 
larger than the nominal uncertainty that Case 
& Bhattacharya claimed (although there is ev- 
idently a wide spread of intrinsic properties of 
SNRs with known distance, see Fig. 0}. If a 
E oc D" relation is to be used to describe the 
relationship between E and D, then arguably 
(see Isobe et al.) a fit to observations that treats 
E and D symmetrically is appropriate (e.g. the 
bisector of least square fits made minimizing 
either log E or log D). 

6. Galactic SNR distribution 

The distribution of SNRs in the Galaxy is of 
interest for many astrophysical studies, partic- 
ularly in relation to their energy input into the 
ISM and for comparison with the distributions 
of possible progenitor populations. Such stud- 
ies are, however, not straightforward, due to 
observational selection effects and the lack of 
reliable distance estimates available for most 
identified remnants. In particular, all SNRs in 
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the anti-centre (i.e. 2nd and 3rd Galactic quad- 
rants) are outside the Solar Circle, at large 
Galactocentric radii. These are regions where 
the background Galactic emission is faint, so 
that low surface brightness remnants are rela- 
tively easy to identify (see SectionfJ}. Without 
taking selection effects into account, the larger 
number of fainter SNRs in the anti-centre leads 
to an apparently broad distribution of Galactic 
SNRs in Galactocentric radius (e.g. see further 
discussion in Green 2004). A simple approach 
is to compare the distribution in Galactic lon- 
gitude of bright remnants, with E\ ghz > 10~ 2() 
W irT 2 Hz 1 sr~' where selection effects are 
not important, with th e expected distribution 
for various models. In iGreenl (120041) this was 
done with 64 bright remnants, and simple 
Monte Carlo models of a Gaussian distribu- 
tion, where the probability distribution varies 
with Galactocentric radius, R, as oc e~ (S/,<r) ~, 
where <x is the Galactocentric scale length. A 
simple least squares comparison 6 of the ob- 
served and model cumulative distributions in- 
dicates that for this simple model a scale length 
of cr ~ 6.5 kpc best matches the observed dis- 
tribution of high brightness SNRs. 

7. Conclusions 

In statistical studies of Galactic SNRs it is cru- 
cial that the selection effects that apply to the 
identification of remnants are taken into ac- 
count. Both faint (i.e. low surface brightness) 
and small angular size remnants are missing 
from current catalogues of Galactic SNRs. A 
consequence of the angular size selection ef- 
fect is that few young but distant remnants have 
yet been identified in the Galaxy. For remnants 
with known distances, the intrinsic range of 
luminosity of Galactic SNRs is large, which 
combined with selection effects, means that the 
L-D relation is of limited use for determining 
distances to individual remnants, or for statis- 
tical studies. 
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6 Not a y 2 comparison, as erroneously stated in 
lGreerJ<2004 . 
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